x1=[4 1 ;2 4;2 3;3 6;4 4];
x2=[9 10;6 8;9 5;8 7;10 8];

s1=5*cov(x1);
s2=5*cov(x2);
u1=mean(x1);
u2=mean(x2);

sw=s1+s2;
sb=(u1-u2)'*(u1-u2);
[v d]=eig(inv(sw)*sb);
v2=inv(sw)*(u1-u2)';

%using LDA
data=[x1;x2];
target=[ones(length(x1),1); zeros(length(x2),1)];
w=LDA(data, target);
